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Abstract. We examine the question of survival of quantum entanglement between 
the bipartite states and multiparticle states like GHZ states under the action of a 
dephasing bath by the application of sequence of tt pulses. We show the great 
i-£h ; advantage of the pulse sequence of Uhrig [ 2007 Phys. Rev. Lett. 98 100504] applied at 

irregular intervals of time, in controlling quantum entanglement. In particular death 
of entanglement could be considerably delayed by pulses. We use quantum optical 
C^i ' techniques to obtain exact results. 
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1. Introduction 

It is well known that the quantum entanglement deteriorates very fast due to 
environmental interactions and one would like to find methods that can save or at 
least slow down the loss of entanglement pQ. It is also now known that quantum 
entanglement can die much faster than the scale over which dephasing occurs [2]. For 
example the coherence of the qubit typically lasts over a time scale of the order of T 2 
whereas the entanglement can exhibit sudden death and thus it is important to extend 
the techniques used for single qubits to bipartite and even multipartite systems. In this 
paper we examine how the pulse techniques which were developed to examine the issue 
of dephasing can help in saving the entanglement. The quantum dynamical decoupling 
[31 SI [5] uses a sequence of control pulses to be used on the system at an interval much less 
than the time-scale of the bath coherence time. In this way, the coupling of the system 
to the bath can be time-reversed and thus canceled. Such non-Markovian approach 
has been successfully applied to two-level systems, harmonic oscillators [6J. A different 
approach was used in [7] where a control pulse was applied to a different transition rather 
than the relevant two-level transition. This technique shows that the control pulse causes 
destructive interference between transition amplitudes at different times which leads to 
inhibition of the spontaneous emission of an excited atom. Similar techniques could 
be useful to suppress the decoherence of a qubit coupled to a thermal bath. Other 
methods for protection against dephasing are known. These include application of fast 
modulations to the bath [8] as well as decoherence free subspaces [9]. The dynamical 
decoupling idea has been implemented in a few recent experiments [10] with excitons in 
semiconductors, with Rydberg atomic qubits, with solid state qubits and with nuclear 
spin qubits. 

More recent developments primarily due to Uhrig [HI [121 H3J EH] go far beyond 
than what has been done earlier on dynamical decoupling. The dynamical decoupling 
schemes use a series of 7r pulses applied at regular interval of times. The pulses reverse 
the evolution given by the Hamiltonian describing the interaction with a dephasing 
environment. This is because under a it pulse the spin operator S z reverses sign. Uhrig 
discovered that 7r pulses applied at irregular intervals of time are much more effective 
in controlling dephasing. The regular pulse sequence and the Uhrig sequence are given 
by 

Ti = -^—^ =Tsin 2 ( * j - ). (1) 
3 n + 1 ' J V 2(n + 1) J V ; 

In this paper we focus on the utility of the sequence of pulses as discovered by 
Uhrig in saving quantum entanglement. Unlike other papers which focus on dephasing 
issues we concentrate on entanglement. This is important as the dynamical behavior 
of entanglement could be quite different than that of dephasing. We calculate the 
concurrence parameter [15] which characterizes the entanglement between the two 
qubits. We show the net time evolution of the concurrence parameter under the action 
of the Uhrig sequence of pulses and compare its evolution with the one given by when 
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the uniform sequence of pulses is applied. We show the great advantage of the Uhrig 
sequence over the uniform sequence in saving entanglement. A very recent experiment 
[16] establishes the advantage of Uhrig's sequence in lengthening the dephasing time of 
a single qubit. The organization of the paper is as follows: In Sec 2 we introduce the 
microscopic model of dephasing and calculate the relevant physical quantities under the 
influence of the control pulses. In Sec 3 we show how the coherent state techniques can be 
used to obtain the dynamical results. In Sec 4 we calculate the dynamics of entanglement 
and present numerical results. In Sec 5 we conclude with possible generalizations of our 
results on entanglement. 



2. Dynamical decay of entanglement under dephasing 



Let us consider two qubits in an entangled state [2] which in general could be a mixed 
state. In terms of the basis states for the two qubits, we choose the initial state as 
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The state ([3]) is positive and normalized iia + b + c + d= l and be > \z\ 2 . This state 
has the structure of a Werner state. For a — d = 0; b = c = \z\ = 1, it represents a 
maximally entangled state. The amount of entanglement in the state is given by the 
concurrence given by 

C = Max{0,(7}, 

C = 2{\p 23 \- y /p^} = 2\z\{l-r), 
\fad 



(4) 



And therefore the state is entangled as long as \z\ is greater than y ad. Now under 
dephasing the diagonal elements a and d do not change. However the coherence in 
the qubit decays as exp[— t/T 2 ] and therefore the entanglement survives as long as 

\z\ exp[— 2t/T-i\ — y/ad > and thus entanglement vanishes if t > — In — 

2 V ad 

We would now examine how the action of pulses can protect the entanglement. 
We would calculate the time over which entanglement can be made to survive. For 
this purpose we need to examine the microscopic model of dephasing. We would make 
the reasonable assumption that each qubit interacts with its own bath. We could then 
examine the dynamics of the individual qubits and then obtain the evolution of the 
concurrence. 
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On a microscopic scale the dephasing can be considered to arise from the interaction 
of the qubit with a bath of oscillators i.e. from the Hamiltonian 

H = h^Uiajai + hS z ^2gi(ai + a\), (5) 

i i 

where the S z is the z component of the spin operator for the qubit and the annihilation 
and creation operators a^, af represent the oscillators of the Bosonic bath. The bath is 
taken to have a broad spectrum. In particular for an Ohmic bath we take spectrum of 
bath as 

J — > X! \9i\ 2 H u ~ u i) = ^olojQ{ijJd — oj), (6) 

i 

where un is the cut off frequency. It essentially determines the correlation time of the 
bath. Such a bath leads to dephasing i.e. the spin polarization decays at the rate T 2 . 
The dynamical decoupling schemes use a series of 7r pulses applied at regular interval of 
times whereas Uhrig applies tt pulses at irregular intervals of time. Such nonuniformly 
spaced pulses are much more effective in controlling dephasing over a time interval 
determined by the cut off frequency and number of pulses. The regular sequence of 
pulses is more effective outside this domain. The pulses reverse the evolution given by 
the interaction part in the Hamiltonian ([5]) since under a tc pulse the spin operator S z 
reverses sign. The regular pulse sequence and the Uhrig sequence are given by equation 
(PO). We need to calculate the dynamical evolution of the off diagonal element of the 
density matrix for the qubit. We work in the interaction picture hence the Hamiltonian 
fl5]) becomes 

H = hS z J2gi(a i e- iUit + a\e^ t )=hS z B(t), (7) 

i 

where B(t) is the bath operator given by 

B(t)=J2g*(a l e-^ t + a\e^ t ), (8) 

i 

It is easy to see that the off diagonal element of the single qubit density matrix a 

is 

a n (t) = Tr B U \U(t)a B a(0)tf(t)\ f), (9) 
where Tr# is over the initial bath density matrix erg and where 

U{t) = Texp{-t [ S z B(t)g!t}. (10) 



This can be simplified to 

vn(t) = ff n(o) Tr B V-(t)<r B vl(t), 



where 



(11) 

a n (0)(Vl(t)V4t)), 



V ± (t)=Te W {T l - f B(r)dr}. (12) 
2 Jo 
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Figure 1. A sequence of it pulses is applied at times Tj. 



Thus we can write 

a n (t) = a n (0)C(t), (13) 

C(t) = (Vl(t)V4t)), 

(14) 

= (Texp{iJ*B{r)dT}). 

So far no approximation has been made. Now we incorporate the effect of pulses in the 
dynamical evolution of the single qubit coherence. Let us apply a sequence of 7r pulses 
at times Tj as shown in figure 1. At each Tj the interaction Hamiltonian changes sign. 
This can be easily incorporated in the dynamics and the result is 

m = (w(t)), 

W(t)=Texp{if*B(T)f(T)dT}, 

(15) 

N-l 

/(t)=E(-W-Ti)^ + i-i), 

j=0 

where the step function 9(t) — 1 if t > 0, = 0, if t < 0. It is especially instructive to 
use coherent state techniques to simplify the expression for W. We do this in the next 
section. 

3. Single Qubit Coherence S(t) 

We now examine the calculation of the function W(t). We note that the bath operator 
B(t) is such that the commutator [B(ti), Bfo)] is a c-number. In such a case it has 
been shown by Glauber [17] that the time ordering can be simplified. It can be shown 
that 

W = exp{z f B(r)drf(r)}exp{-1- f d n P dr 2 [S(n), S(r 2 )]/(ri)/(r 2 )}. (16) 
Jo z Jo Jo 

Since B is a Hermitian operator, the last exponential is just a c-number phase factor 

and hence 

W = exp(i$(t)) exp{i f BMdrfM}, (17) 

Jo 

= exp(i$(t))IIj expiifjOj + if* a]}, (18) 
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fj=9j f e-^ T f{r)dr. (19) 

J 

On using the Baker-Hausdorff identity (|T8|) can be further simplified to 

W = exp(^(t))n,exp{^;4}exp{^a :; }exp{-^|/ J | 2 }, (20) 

The thermal expectation value of W can be easily obtained using for example the P- 
representation for the thermal density matrix [18] 
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Pthj = / exp{ }|a)(a|d a 



(21) 



Thus 

W = exp(^(t))n, exp{-i|/ i | 2 } x — / exp{-^L} exp{z/*a* + if ja }cPa, (22) 

2 Tinj J rij 

which on simplification reduces to 

W = exp(i$(t))n i exp{-(n i + \)\f 3 \ 2 }, (23) 

On using the form of fj and on introducing the spectral density of the bath oscillators, 
the expression (123]) becomes 

W = exp(i<&(t))exp{- J < W )JMK W ) + i]|/( W )| 2 }, (24) 

where now 

f(u)= fe-^F{r)dr. (25) 
Jo 

The function / can be simplified using the explicit form of F(r): 

N 

f(u) = -i[l + (-l) N+l e- iuJt + 2j2(-^ J e-~ iujT '}. (26) 

j'=i 

The result ( 1231 is equivalent to equation (jS|) of Uhrig. We also note that results like 
(24) appear in the earlier literature (8] dealing especially with nonmarkovian master 
equations. 



4. Saving Entanglement: Numerical Results 

Since we work under the assumption that each qubit interacts with its own bath, the 
time dependent matrix elements of the density matrix in the basis ([2]) can be obtained 
by noting that the diagonal elements do not evolve under dephasing. The off diagonal 
element is given by 

p 23 (t)=P23(0)|C(t)| 2 , (27) 
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Figure 2. Signal vs time for n = 10 at T — 0. Solid lines for the optimized 
sequence, dotdashed lines for the equidistant sequence. From bottom to top the curves 
correspond to a = 0.25,0.1,0.01,0.001. 




Figure 3. Signal vs time for n — 50 at T — 0. Solid lines for the optimized 
sequence, dotdashed lines for the equidistant sequence. From bottom to top the curves 
correspond to a = 0.25,0.1,0.01,0.001. 



where ((t) is defined by equation f|T3|) and its explicit form is given by equation fl23|) . 
Thus 

|C(t)| 2 = S(t) = exp{-2 J dwJ(u)(n(u) + ^)\f(u)\ 2 }, (28) 
It can then be shown that the time dependence of the concurrence is given by 
C(t) = Max{0,C(t)}, 
C(t) = 2\z\{S(t)-r}, 



ad 



\z\ 

We next discuss the dynamical behavior of the entanglement. The function f(u) 
has been evaluated by Uhrig. For the pulses applied at regular intervals and for n even 
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we have 

\f{oo)\ 2 = 4tan 2 [wt/(2n + 2)]cos 2 (wt/2)/w 2 V n even, (30) 

whereas for the Uhrig's pulse sequence 

\f(u)\ 2 ^lQ(n + l) 2 J 2 +1 (co/2)/u J 2 , (31) 

where J n is the Bessel function. The function S(t) is shown in figures 2 and 3 for n = 10 
and n = 50. The parameter u/^ 1 is a measure of the bath correlation time. These figures 
show that the entanglement lives much longer for Uhrig sequence of pulses applied at 
nonuniform intervals of time provided that uont < In. Thus entanglement can be made 
to live over times which could be several orders longer than the coherence time of the 
bath. 



5. Conclusions 



In conclusion we have considered how the effects of dephasing on the destruction of 
entanglement can be considerably slowed on by applying the sequence of n pulses applied 
at time intervals given by Uhrig. We demonstrated this explicitly for the case of a mixed 
entangled state of two qubits. The sequence given by Uhrig is far better in controlling 
the death of entanglement compared to the sequence applied at regular intervals of time. 
These conclusions also apply to the multiparticle entangled state like the GHZ state 

l*> = ^(lt"-t>-U-"±», (32) 

whose entanglement under dephasing would decay as the density matrix at time t would 
be 

= t • • • tXt • • • 1 1 + • • • iXi • • • i I 

1 tN 
-^x P {- — K\t---t)(i---i\+cc). 

Under the application of ir pulses, the prefactor exp{ — — -} would be replaced by 

T 2 

(S(t)) N / 2 . Since S(t) can be made close to unity for times even orders of the correlation 
time of the bath, the entanglement of the multiparticle GHZ state would survive 
over a long time. Note further that Uhrig's work has been generalized to arbitrary 
relaxations [12J. Clearly these generalizations should be applicable to the considerations 
of entanglement. In particular we hope to examine the protection of Werner state against 
different models of environment. 

Finally we note that our ongoing work also suggests how other methods like 
photonic crystal environment can be used to save entanglement. 
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